On the low-energy density of states in disordered d-wave superconductor 
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We study the low-energy density of states of Dirac fermions in disordered d-wave superconductor. 
At zero energy, a finite density of states is obtained via tlie mechanism of dynamical mass generation 
in an effective (l-l-l)-dimensional relativistic field theory. 
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The low-energy excitations of d-wave cuprate super- 
conductor are relativistic massless Dirac fermions. If the 
superconductor is clean, the density of states (DOS) van- 
ishes as N{uj) ~ in the vicinity of Fermi surface. 
This behavior is believed to be substantially influenced 
by weak static disorder. However, despite extensive ex- 
perimental and theoretical investigations, the low-energy 
DOS of disordered d-wave superconductor is still a sub- 
ject under fierce debatei^^>^>&>2>a>a>io^a2^^^. 

In the conventional self-consistent T-matrix approxi- 
mation (SCTMA), the single-particle density of states 
(DOS) N{lo) at Fermi surface is found to be finite as 
the result of impurity scattering^ii^iiii^i^. The finiteness 
of zero-energy DOS iV(0) received support from early 
numerical calculationsliSii^. However, the prediction of 
finite DOS at Fermi surface and even the SCTMA treat- 
ment were criticized by Nersesyan et al.^. Describing 
the static disorder in terms of random abelian (or non- 
abelian) gauge fields and then averaging the random vari- 
ables, they mapped the two-dimensional disordered Dirac 
fermion system to an effective (l-l-l)-dimensional rela- 
tivistic field theory. At zero energy a; = 0, the effec- 
tive theory possesses a continuous chiral symmetry, with 
iV(0) being the relevant order parameter. Nersesyan et 
al. argued that N{0) must be zero because the continuous 
symmetry can not be spontaneously broken, as required 
by the Coleman-Mermin- Wagner (CMW) theorem. They 
then applied the non-perturbative bosonization method 
to treat the effective theory and found that N{u}) has a 
power-law dependence on energy and finally vanishes at 
w = 0. 

After the work of Nersesyan et al., a series of novel 
methods, mostly non-perturbative, have been put for- 
ward to calculate the low-energy DOS ^^i^^'^'^ . Surpris- 
ingly, the resulted theoretical predictions are quite dif- 
ferent. The DOS at Fermi surface is found to vanish or 
diverge, depending on the disorder model studied and 
methods used. These seemingly conflicting results, how- 
ever, can be unified by introducing a novel diffusion 
mode at momentum (tTjTt), which arises from a global 
nesting symmetry^^. In addition, a duality relation was 
proposed^^ to exist between the Gaussian limit of weak 
disorder and the unitary limit of dilute disorder, thus 
connecting the vanishing A^(0) to the diverging one. 

Though making significant progress in enriching our 
understanding on random Dirac fermions, the theoret- 
ical predictions of vanishing and diverging zero-energy 



DOS are unlikely to agree with experiments in cuprate 
superconductors^. On the contrary, there are experimen- 
tal facts that imply the finiteness of A^(0)^. Specifically, 
the heat transport measurements^^ provide rather strong 
evidence for the existence of a universal thermal conduc- 
tivity, which was deduced from a finite A^(0)^. 

In this paper, we revisit the issue of low-energy DOS 
of disordered Dirac fermions. For weak static disorder, 
the frequencies are not mixed during the process of im- 
purity scattering. Then the frequency ujn can be consid- 
ered as an external parameter. In the effective (l+I)- 
dimensional relativistic model in terms of fermionic field 
77, ll}„ plays the role of a bare mass. It is renormalized by 
random impurity effects to become a physical frequency. 
Due to the mechanism of dynamical mass generation, the 
DOS N{u!) has a nonzero magnitude even at w = 0. 

In the two-component Nambu representation <f>j^ = 

(cj^l, c_ki), the Hamiltonian of two-dimensional d-wave 
superconductor can be written as 



k 



AkT2]$k, 



(1) 



where t, are the Pauli matrices. The electron disper- 
sion is Ek = —2t{coskxa -\- cos kyu) — fi and the d-wave 
energy gap is Ak = Ao(cosfc xd cos/byfl) which has 
four nodes at the Fermi level. Making a 7r/2 rotation 
in the Nambu space as $k e'^'^^^'^^k and then lin- 
earizing the spectrum in the vicinity of nodes, we have 

_Ek — ±1 



,2 1,2 



■jfr^i -r V2k2, where ki is perpendicular to the 

Fermi surface and k2 is parallel to the Fermi surface. We 
also defined the Fermi velocity "Vf — de^^/dk — Vfki with 
Vf ^ 2\/2ta, and the gap velocity V2 = 9Ak/9k = V2i^2 

with V2 = Aoa/\/2. 

Following Nresesyan et alr^, we use a continuum de- 
scription in terms of the two-component spinor field ip 



Hi 



(2) 



Hu2) = -Wl9l(2)Tl - iV2d2(l)T2, (2) 



with node index being j = 1,1,2,2. The next step 
is to rescale the spatial coordinates and spinor field as 
xi VfXi, X2 V2X2, and ijjj {vfV2)~^^'^4'j ■ For the 
sake of simplicity, we consider only the forward scatter- 
ing process, so the random potential couples to the Dirac 
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fermion as a random abelian gauge field-lS. For the dis- 
ordered Dirac fermion with frequency a;„, the generating 
functional is written in the form 



Z = 



(3) 



with two-dimensional Euclidean action 

S = \ d^xtp^ [-iTi {di + iA) - iTa^a - iujn]^p. (4) 



The probability distribution of random gauge field A{x) 
has the form 



P[A] ^ J VAexp j (f 



A^jx) 
' 2.9 



(5) 



Then the random gauge field can be averaged using the 
replica trick, leading to the following effective action 

5 = 1 d'xmrd, - u;„)v - ^{fil'vn (6) 

where the new fermion fields are defined as fj = — z-^^ 
and rj = ijj. This action has already been written in the 
(l+l)-dimensional Minkowski space with metric g^^ = 
(1. — 1). There is a summation over all r replicas, though 
the summing notation is hidden. The gamma matrices 
are given by 7° = T2 and 7^ = ti. 

The action (6) is identified as the (H-l)-dimensional 
massive Thirring model, and frequency ujn as the mass 
of fermionic field 77. Within this model, the disordered 
Dirac fermion DOS is 



N{uj) - Re[Tr(^(a;)77(a;))]| 



iiOn — >uj-\-i6 • 



(7) 



At zero-energy = 0, the action (6) is invariant under 
the continuous chiral transformation 



V 



rje 



(8) 



where 7^ = ra. A mass term ^ {fjr]) seems to break this 
symmetry. It was claimed^*^ that the zero-energy DOS 
A^(0) must vanish because continuous symmetry can not 
be broken in (l+l)-dimensional. However, the CMW 
theorem does not apply to the present case since the 
associated symmetry group is non-compac t^^'^^i^^ . The 
massless fermion rj can acquire a finite mass through self- 
interactionSS. Dynamical mass generation for fermion 7/ 
should have an important effect on the low-energy DOS 
of the massless nodal fermion. 

Now we calculate the dynamical mass by solving the 
self-consistent Dyson- Schwinger integral equation. In the 
effective (l-l-l)-dimensional model (6), the free fermion 
propagator is 



Go(fc) 



7 • fc - 



(9) 



The bare mass a;„ is corrected by the disorder effects to 
become a new quantity a3„, which may be called a phys- 
ical mass. The associated complete fermion propagator 



can be written as 



G(fc) = 



7 • - u;„ 

The Dyson-Schwinger equation is given by 



(10) 



G-\p) = G,\p)~Tt 



{2n) 



Gik)T{k,p). (11) 



To the one-loop order, the vertex function T(k,p) is sim- 
ply the bare vertex —iX — —igN, where N is the number 
of nodes (flavor of Dirac fermion). Substituting propa- 
gators (9) and (10) into the above equation, we have 



— UJn + A 

At zero frequency aj„ 
ing solution 

',0 — 



_fk 
(2^ 



rTr 



7 • fc 



(12) 



0, the equation has the foUow- 



w„ = 7 = Aoe" 



(13) 



For finite but small frequency aj„, w„ is the frequency- 
dependent, physical mass. It is interesting to point out 
that this kind of nonlinear equation has been encountered 
in the QCD theory^. Through the strong interaction 
mediated by non-abelian gauge fields, the small bass mass 
of quarks is significantly enhanced via the mechanism 
of dynamical (approximate) chiral symmetry breaking^^. 
After angular integration and Wich rotation, the integral 
equation becomes 



UJ„ 



1 



(27r)2 fc|. 



~ A Ao 

- LL>n — in -zz— ■ 



(14) 



where Aq acting as the ultra-violet cutoff. Suppose that 
LOn <C Wn, then the solution is 



7 + 



A 



(15) 



which has the correct limit 5„(a;„ — > 0) = 7. When 
the energy ojn > 7, the quantity 7 is relatively unim- 
portant and can be safely neglected. As the frequency 
is decreasing, it becomes comparable with ujn and finally 
dominates at a;„ ^ 0. 

To calculate the low-energy DOS of fermions, we may 
return to the (2-|-l)-dimensional Matsubara formalism. 
The DOS is given by 



N{uj) = 

TT 



(2^)2 

Using the retarded Green function 
G'ret(k,w) 



Tr[ImG„t(k,w)]. (16) 



— T • k + 



(17) 
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we arrive at the following low-energy DOS 



N{uj) 



7 



TTVfV2 VfV2 



X' 



(18) 



At zero energy limit, it has a finite value iV(0) = 
j/TTVfV2- This result can be understood as follows. As- 
sume that the renormalized frequency w contains all the 
effects of random impurity potentials. Then the disor- 
dered Dirac fermion with frequency lo is equivalent to 
the clean Dirac fermion with frequency uj. It is known 
that the DOS of clean Dirac fermion is^ 



TTVfV2 



(19) 



Replacing bare frequency w with dynamical lo then leads 
to the DOS of disordered Dirac fermion (18). 



In conclusion, while starting from the same Hamilto- 
nian and effective (l-l-l)-dimensional field theory as those 
of Nersesyan et aL^°, we arrive at a different expression 
of low-energy DOS, which remains finite at zero-energy. 
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